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A STUDY OF THE STABILITY OF THE LAMINAR BOUNDARY LAYER 
AS AFFECTED BY CHANGES IN TEE BOUNDARY -LAYER 
THICKNESS IN REGIONS OF PRESSURE GRADIENT 
AND FLOW THROUGH THE SURFACE 
By Neal Tetervin and David A. Levine 


SUMMARY 


The Schlichting method for the computation of the laminar boundary 
layer states that on an impervious surface the shape of the velocity pro- 
file is determined by the local effective pressure gradient) this gra- 
dient is directly proportional to the product of the actual pressure gra- 
dient and the square of the boundary-layer thickness. It is apparent, 
therefore, that in a region of falling pressure the Schlichting method 
predicts that an increase in boundary -layer thickness increases the 
effective pressure gradient and thus results in a more convex velocity 
profile. Because the increase in convexity is known to imply an increase 
in the critical boundary-layer Reynolds number, it would appear that an 
increase in boundary-layer thickness could increase the local critical 
Reynolds number more than the local boundary-layer Reynolds number. 

In order to investigate this possibility, computations have been 
made by combining the Schlichting method with the Lin method for the 
calculation of the critical Reynolds number of a velocity profile. The 
computations indicate that in a region of falling pressure on an imper- 
vious surface an increase in boundary -layer thickness can cause the 
velocity profile shape to be changed enough by the increase in effective 
pressure gradient so that the ratio of the local critical Reynolds num- 
ber to the local boundary -layer Reynolds number is increased. It thus 
appears that the local stability of the boundary layer can be increased 
by a local increase in boundary -layer thickness. The computations also 
indicate that similar effects occur when there is flow through the sur- 
face) in this case the results depend on the effective flow through the 
surface as well as on the effective pressure gradient. 

These calculations suggest that an increase in boundary-layer thick- 
ness can decrease the disturbing effect of roughness particles without a 
decrease in stability. This conclusion is implied by the result that an 
increase in boundary-layer thickness reduces the velocity at a fixed 
distance from the surface more than the change in velocity profile 
increases the velocity. One method of increasing the boundary-layer 
thickness, namely, blowing near the stagnation point, has been investi- 
gated theoretically and seems to have limited potentialities because the 
blowing produces a significant increase in boundary-layer thickness only 
over the foremost portion of the airfoil. 
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INTRODUCTION 

A number of years ago investigators noticed that a falling pres- 
sure in the direction of flow was able to increase the extent of laminar 
flow on a surface and thereby decrease its friction drag (references 1, 

2, and 3)* This observation was made the basis for the design of air- 
foils that are able to have much smaller drag coefficients than previ- 
ously known airfoils. In order to -- obtain the small drag coefficients, 
however, it is necessary that— the airfoil surface be smooth and free of 
waves. The requirement of “a smooth and wave -free surface has, to the 
present, prevented the consistent attainment in flight of wing drag 
coefficients as small as those obtained in wind-tunnel tests. 

The experimental finding that under proper conditions a falling 
pressure in the direction of flow causes increases in the extent of 
laminar flow was given a theoretical basis by the work of reference 4. 

The theoretical work showed that a pressure drop increases the convexity 
of velocity profiles and thereby Increases their critical Reynolds num- 
bers. The critical Reynolds number is the Reynolds number below which 
transition to turbulent flow cannot - be caused by the growth of small 
departures from the mean velocity inside the boundary layer as the fluid 
moves downstream. To form the Reynolds number the velocity at the outer 
edge of the boundary layer is used; the length is a boundary-layer 
thickness . 

Another important theoretical finding (reference 5) was that a flow 
of fluid -into the-surface increases the convexity of the velocity profile 
and thereby increases its - critical Reynolds number; an outflow was found 
to decrease the - critical Reynolds number. Because the effects on the 
critical Reynolds number are large for very small flow velocities through 
the surface, the method of maintaining laminar flow by sucking fluid into 
the airfoil through a porous surface appeared to have practical value. 

In order to test the method, the work of reference 6 was done. Full-chord 

laminar flow was observed up to a Reynolds number of 20 x 10^ and it was 
concluded that larger Reynolds numbers would" not limit the extent of lami- 
nar flow if the surface were sufficiently smooth and free of waviness. 

Another method for maintaining laminar" flow was based on the experi- 
mental observation that transition often seems to occur at a fairly defi- 
nite value of the boundary-layer Reynolds number. Examples of the appli- 
cation of this observation are references 7> 8, and 9 in which a number 
of slots were placed along the airfoil surface and air drawn into the 
interior in order to decrease the boundary -layer Reynolds number. 

In tests of each of the three methods, namely, pressure drop, suc- 
tion through a porous surface, and limitation of magnitude of boundary- 
layer Reynolds number, it was found that, in' order to avoid early transi- 
tion to turbulent flow, the surface has to be free of noticeable roughness 
particles and other departures from smoothness. The work of reference 10 
indicated that the roughness Reynolds number, formed from the height of- 
the roughness particle and the velocity at - a“ : distance from the surface 
equal to the roughness height with the particle absent; furnishes a 
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> measure of the disturbing effect of a particle. For a small effect, a 

small roughness Reynolds number is necessary. All three methods for 
keeping the boundary layer laminar, however, increase rather than 

v decrease the roughness Reynolds number of a particle. 

In the Schlichting method (reference 11), the shape of the velo- 
city profile on an impervious surface depends only on the kinematic 
viscosity, on the local gradient of velocity at the outer edge of the 
boundary layer, and on the local boundary-layer thickness. The combina- 
tion of variables is such that an increase in boundary-layer thickness 
increases the effect of velocity gradient; therefore, in a region of 
falling pressure, an increase in boundary- layer thickness results in a 
more convex velocity profile. Because the increase in convexity carries 
with it an increase in the critical boundary-layer Reynolds number, the 
possibility arises that an increase in boundary-layer thickness may 
increase the local critical Reynolds number more than the local boundary- 
layer Reynolds number is increased. by the increase in thickness. There 
is thus the possibility that an increase in boundary- layer thickness can 
increase the stability of the boundary layer. 

1 The purpose of the present work is therefore to investigate theo- 

retically the effect of an increase in boundary- layer thickness on the 
stability of the boundary layer. The analysis also includes an investi- 
gation of the effect of an increase in boundary- layer thickness on the 
roughness Reynolds number and of the ability of blowing at the stagna- 
tion point to increase the boundary-layer thickness without making the 
boundary layer ms table. 


SYMBOLS 


a,b coefficients in expression for t|t 



c Dp 


drag equivalent of power required for boundary- layer control 


sum of wake drag coefficient and Cpp 



c chord of airfoil 


c 0 ,c-[_,C2 coefficients in expression for g 
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efficiency of boundary- layer control system 
efficiency of aircraft - propuls ive system 
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*1 - - 

“oV V 


P-L = 1 - e -r l 


F 2 = 1 - e _1 l - sin g tj 


(0 < n < 3) 


= -e -T l 


function of" y 


(n > 3) 


s = zr- = c Q -t-c-jK + CpK 2 
8l 

S = 2 f - k^2 + ■ k L 


height of roughness particle 
total pressure 


v^c *i 2 . 

2 dU 2 — 2k 
3x dx ' 


velocity-profile shape parameter 
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static pressure 

free -stream static pressure 

power required for boundary- layer control 


volume flow per unit time 


critical Reynolds number, value of Rg at which, a small 
disturbance is neither damped nor amplified 


critical Reynolds number, value of Rg* at which a small 
disturbance is neither damped nor amplified 


roughness Reynolds number 


airfoil Reynolds number 


velocity at outer edge of boundary layer 
velocity in direction of x 


u = iL 


free -stream velocity 
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Ujj velocity at y = h with, roughness particle absent 

v w velocity through surface, positive outward , v 



u 0 

v velocity normal to surface, positive outward 

V velocity in direction of y-axis, stagnation-point flow 



x distance along surface measured from stagnation point 



c 


Ax increment in x 

y distance normal to surface, positive outward 



c 


Z = 0 2 R c 

a angle of attack 

8^ measure of boundary- layer thickness 
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momentum thickness 


Q2 3TT R 


k = §1 £1 = l£L/_l_ z £ 

V d5c R c \U 2 dx/ 


k l = - 


- — R e = 
u 


x = h-^ 

v dx 


X 1 = ~ 


Vi 


viscosity of fluid 

kinematic viscosity of fluid (TT/p") 

free-stream density 

density of fluid 

surface shearing stress 


function of y 


nondimens ional stream function 


stream function 


Subscripts: 


stagnation point 


at surface 
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00 

far behind wing 



X 

at station x 



x+Ax 

at' station x + Ax 

— 



Barred quantities are dimensional. 


ANALYSIS 


Effect of an increase in boundary- layer thickness on the stability 
of the laminar “boundary layer . - The analysis makes use of a modifica- 
tion of Schlichtlng’s method for the calculation of laminar -boundary- 
layer velocity profiles (reference 11). The modification consists in 
using the boundary-layer momentum thickness 6 as the reference length 
instead of the length used by Schlichting. The use of G instead 

of — 5p eliminates two auxiliary parameters, X and and makes the 

extension of the Schlichting method to cases of uniformly distributed 
blowing relatively simple. The modified Schlichting method is described 
in appendix A. 


From the definition of the critical Reynolds number (reference 12) 
it follows that the boundary layer is unstable when the ratio of the 
local critical Reynolds number- to the local boundary- layer Reynolds 


number 


— — is less than unity, neutrally stable when = 1, and 

Ra Re 


stable when 


R, 


6 


the numerical value of- the ratio 


>1. In the present work the assumption is made that 

R 0c 


Re 


at a distance x from the 


stagnation point - is a measure of the stability of the boundary layer 
at that .value of x. It is assumed that the” boundary layer becomes 

R 0 

more stable as the ratio — — increases. (in this analysis it is 

Re 

more convenient to use Rg and R q than Rg* and Rg* c .) 


A conclusion of reference 13 is that- the critical Reynolds number 
Rq c depends only on the shape of the velocity profile. In the 

Schlichting method all velocity profiles are a single -parameter family. 
Therefore to each value of- K, the velocity-profile parameter, there 
corresponds a value of Rg c ; that is, 

Rg c = Rg c (K) 


£ 


V 


( 1 ) 
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The function. Rg c (K) was calculated "by Lin’s rapid method (reference 12) 
and is shown in figure 1. (See also reference 6. ) 

The possibility that an increase in 9 can increase K enough 
to cause a larger increase in R q c than in Rg is now investigated. 

The investigation is made for the general case v^ f 0; v w = 0 is 
then a special case. 


In order to determine the effect on 


Re 


at a specific value of x 


- — dU _ 

of an increase in 6 at that value of x with v. U, — , and v™ 

dx w 

R 0c 

Re 


fixed ^R c , U, — , and v w fixed^, it is convenient to obtain 


dx 

a function of the parameters K and 


as 


v w 2r c 

2 M ' 

dx 


The shape parameter K is 


a function of the two quantities k and kp defined by: 


e£ dU 
v dx 


and 


ki = 



These definitions can also be written as: 


_ R e 2 1 dU 

R c U 2 dx 


( 2 ) 


and 


ki = - — R 0 

u 


( 3 ) 
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When the definitions for k and kq are substituted into the 
following equation connecting K, k } and :kq (equation (A8 ) , 
appendix A): . 


g2(K + 1) - fkq - k = 0 


( 4 ) 


the resultrts 


JL _ f Zw r g 2 (k + 1 ) = 0 


R c U 2 dx U 


( 5 ) 


The quantities f and g where 


f = 


T w 9 

Tiu 


and 


9 

8 = — 
&1 


are functions of K alone. Equation (5) is solved for R§; the result- 
is 


^ R e = Jf 
U 0 


1 ± l/l + 


2 g 2 (K + 1 ) 
Jf 2 


( 6 ) 


where 


J = 


v w R c 
2 

dx 


* 1 2 


-k 2 
K 1 


2 ^ 2k 


dx 
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It is noted that the parameter J is independent of Rg and therefore 
does not change when 9 at a fixed value of x is changed. 

The sign of the radical in equation (6) is determined hy the fact 
that Rg is always positive. The sign is given in the following table: 


J 

K + 1 

Sign 

— > 0 
u 

>0 

>0 

+ 

>0 

<0 

± 

<0 

<0 

- 

v w 

— < 0 
u 

>0 

>0 

- 

<0 

<0 

+ 

<0 

^0 

+ 


By combining equations (l) and (6) the following result is obtained 




Ra c (K) 


— 

u 


jf 


1 ±\/i + 


2g 2 (K + 1) 

Jf 2 J 


(v w ± 0) 


( 7 ) 


R, 


Equation (7) fixes — — when J, and K are known. The variation 

R 0 U 


R 


e. 


V W Rn 


Of 


with K 


for constant values of J 


is shown in figure 2. 
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(The curves of fig. 2 have been numbered; the value of J to which the 
numbers correspond may be obtained from table I. ) 

Usually the airfoil does hot have a porous surface; in this case 
v w = 0 and equation (7) is replaced by 



Equation (8) 
solving for 

as v w — >0 

figure 3. 


can be obtained either by putting v w = 0 in equation (5), 
Rq, and combining with equation (l), or by taking the limit 


in equation (j). 


The function 


He, 

R a 


u 2 R c 

dU 

dx 


is shown in 


The quantity under the square-root sign in equation (6) must be 
positive; thus, 


1 + 


2g 2 (K + 1) 


jr 


> 0 


When the equality sign is used, the result is the curve defined by the 
equation 




Re c (K) 

2x 2 (K + 1) 
f 


(9) 


and shown by dashed lines in figures 2(a) and 2(b). 


A curve of 




for a fixed value of J cannot cross the dashed curve. 
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Although the relation between K, k, and kq (equation (4) an d 
fig. 4) is used in obtaining equations (7) and (8), not all combinations 
of values are allowable because some correspond to a physically unreali- 
zable behavior of the boundary layer. Thus, in order to avoid the con- 
tradiction that, at a fixed value of kq and with k negative < 0^, 

an increase in k (a decrease in the magnitude of — for — < cA can 

^ dx dx y 

result in separation, only the portion of the curves for K, k, and kq 

in figure 4(a) for which (— \ >0 is used. The boundary curve, 

W k q 

shown as a dashed line in figure 4(a), passes through the points at which 
changes from 00 to -<». Its equation is obtained by applying the 

condition 


>(f£) >0 


Wkq 


to the expression. 



(10) 


obtained from equation (4). The subscript on the partial derivative 
denotes the quantity held constant during the differentiation. From 


equation (10), the fact that 
the result is obtained that: 


— > 0, and the 
dK 


inequality for 



dg2(K + 1) 


kl< 


dK 


df 

dK 



NACA TN 2752 


14- 


On the boundary curve the equality sign applies. When the condition 


dg2(K + 1) 



dK 


is used with equation (4), the equation of the boundary curve is found 
to be 


dg 2 (K + 1) 

k = g2(K + l) - f 1| (11) 

dK 

All physically meaningful values of k and 'K lie to the right of the 
dashed boundary curve (fig. 4(a)). For kp = 0 (that is, v w = 0) only 
values of K greater than -1.913 are allowable. 

R 0c 

The curves of in figure 2 thus satisfy the conditions that 

— ^ R 0 
U 9 


R 0 > 0 


1 + 


gg^K + 1) > Q 
Jf2 


and 
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The curve of 



(fig. 3) satisfies the conditions that: 


R 0 >0 



Before the effect of an increase in 9 on the ratio 



can he 


found from figures 2 and 3, the effect of the increase in G on K 

must he found. The sign of the derivative |^— is thus required. 

dR 0 

Only 9 varies; therefore, Rg varies, hut R c , U, — , and v w do 

dx 

not. From equation (k) K = K(k,k L ); therefore. 


3k = /3K\ 3k_ + /3 k\ ^1 

dR e \Skj kl SR e 3 r 0 


where the subscripts denote the quantity held constant. Now from 
equation (2) 


% 


3k 

3R 0 


2k 


R £ 


3k} 

3Rfl 


= ki 


and from equation (3J 
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Therefore, 


R 


SK 
9 Sr 


e 




( 12 ) 


From equation (4) there follows 


and 




1 

dg 2 (K +• 1) 

- M kl 

dK 

dK kl 

f 


dg2(K + 1) 

-a?*, 


dK dJC 


When these expressions for 



into equation (12)^ the result is 


and 



are substituted 


% 


= 2g g (K + 1) - k]f 


SK 


SR0 dg 2 (K + 1) 
dK 


( 13 ) 


- ki 


df - 
dK 


The curves in figure 2 are identified. by the value of the param- 
eter J. (See table I.) Therefore, in order to find the sign of 
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for any curve in figure 2^ it is necessary to introduce the parameter J 


into equation (13) for Rg 7- — . Equation (6) is used to obtain: 


SR 


kl = > _K R 0 = -jf 
U 


1 + 1 + 


2g 2 (K f 1) 


Jf c 


Equation (13) thus becomes: 


R e 


Sk 

Sr 6 


2g2(K + 1) + f 2 J 

l 1 1/1 + 2g2(K + 11 

1/ Jf 2 

dg2(K + 1) , df & t 

!— — 

1 ± 1/1 + 2e2(K + 1) 

dK dK 

1/ Jf 2 


(1A-) 


The quant ity Rg 


SK 

SRq 


is 


of J in figure 5(a) for v w 

_ B Sk 

curves of Rg 


SR, 


against K 


e 


plotted against K for constant values 
> 0 and in figure 5(b) for v w < 0. The 
also satisfy the conditions that 


R 9 > 0 


and 


1 + 


2S 2 (K + 1) > 

Jf 2 



> 


0 
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When. 

then, 


1 + 

Re 


2g 2 (K + 1) 
Jf 2 

8k 


0, the numerator of equation (l4) is zero; 


The curve of R fl 
equation (l4) with J = 0, 
the conditions that R g > 0 


against K for_ v w = 
is shown in figure 6. 



0, calculated from 
This curve satisfies 


Wow consider figures 3 and 8 which are for v w = 0. It is apparent 
from figure 6 that when v w = 0 an increase in Rg decreases K if_ 

K < -1, does not change K if K is initially equal to -1, and 
increases K if K > -1. Therefore an increase in F, and thus in Rg, 


causes 



to increase if — K > -0. 925 


(fig. 3). 


Consequently, an 


increase in boundary- layer thickness increases rather than decreases 
the stability of the boundary layer if K > -O.925. When K < -O.925 
a small increase in boundary- layer thickness decreases the stability 
of the boundary layer; conversely, when K < -0. 925 a small decrease 
in thickness increases the stability. 


For K > -1 and v w = 0, the pressure falls along the surface in 
the direction of flow; when K < -1 and v w = 0, the pressure rises in 
the direction of flow (equation {k)). The values K = -1 and v w = 0 
correspond to a flow with zero pressure gradient. 

Figures 2(a) and 5(a) (for v„ > 0) indicate that when Rg < 0 

dRg 

(curves 6 to 13 inclusive) an increase in Rg decreases the stability. 

When Rg >0 (curves 1 to 5 inclusive), the stability may either 

be increased or decreased by an increase in Rg depending on the values 
of K and J and the magnitude of the increase in Rg. 


For v, 


w 


when R, 


8k 

Rg decreases 


< 0 the conclusion is similar to that for v w > 0; namely, 

< 0 (curves 9 to 11, figs. 2(b) and 5 (b)) an increase in 

the stability. When Rg — 

6Rg 


> 0 (curves 1 to 7 inclusive. 
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figs. 2(b) and 5(t>)) "the stability may either be increased or decreased 
by an increase in Kg, depending on the values of K and J and the 
magnitude of the increase in Rg. For curve 8 the stability is decreased 
by an increase in Rg. 

The computations thus indicate that under suitable conditions an 
increase in boundary-layer thickness can increase the stability of the 
laminar boundary layer. 

Effect of an increase in boundary-layer thickness on the roughness 
Reynolds number R^. - Because the disturbing effect of a roughness 
particle is decreased when its Reynolds number R^ is decreased suffi- 
ciently (reference 10), it is desirable to examine the effect on R^ of 

an increase in 0, that is, an increase in Rg with R c , U, and v w 

fixed. Although an increase in boundary-layer thickness places the rough- 
ness particle at a smaller fraction of the boundary -layer thickness from 
the surface and consequently in a region of smaller velocity, the 
increase in thickness also changes the velocity profile and therefore 
the velocity at a fixed distance from the surface. Whether the increase 
in velocity associated with the change in velocity profile can be 
greater than the decrease in velocity caused by the increase in the 
boundary- layer thickness is now investigated. 

The definition of R^ is 



(15) 


Although the assumption of reference 10 that 



becomes less accurate as the curvature of the velocity profile increases, 
it can be shown that this assumption overestimates the change in veloc- 
ity caused by the change in velocity profile. The assumption is there- 
fore retained in the present analysis. Therefore, equation (15) becomes 



v 
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"but from reference 11 



rf(K) 


Therefore , 


R h . ^ fM 

9 V 


or 


R h = 


h 2 R c 2 U 2 

~G 


f(K) 


(16) 


Because f(K) changes with Rg, it - is not obvious that an increase 
in Rg always decreases R^. The change in R^ caused by a change in 
Rg alone, really ' U alone, is given by: 


6% o P o 

- — = h 2 R c 2 b 2 
aRg 



_L ctf(K) 

R 0 dK 5R e 


or 


^ = V-i + i " Be &) 

aRg Rg\ £ dK 9 aRg / 


When equation (13) for 


Rc 


ax 

aR f 


is used, the_result is : 
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Re 5% 

®h SRg 


dg^(K + 1) 

dK + 2g g (K + 1) 

df/dK f 

dg2(K + 1) 


df/dK 


(17) 


The partial derivative 
for four values of k-^. 


Re ^h 

Rh SRg 


is shown in figure 7 as a function of K 


dR h 

The partial derivative — — is always negative 

^e 


because the numerator of equation ( 17 ) is negative for all values of K 
and the denominator is positive for all values of K. T ha t the denomi- 


nator is positive follows from the condition that 



this 


requires that: 


*1< 


dg2(K + 1) 

dK 

df 

dK 


Therefore, an increase in Rg always decreases R^j that is, an 

increase in boundary-layer thickness reduces the velocity at a fixed 
distance from the surface more than the change in velocity profile 
increases the velocity. Consequently, a sufficient increase in boundary- 
layer thickness always decreases' the disturbing effect of roughness 
particles. 


Effect of an increase in boundary-layer thickness on the stability 
of the boundary layer df two typical airfoils . - The computations thus 
indicate that a sufficient increase in boundary-layer thickness can 
decrease the disturbing effect of roughness particles and that under 
suitable conditions an increase in boundary-layer thickness can increase 
the stability of the laminar boundary layer. An equivalent statement of 


these two results is that 



dRg 


be positive. 




is always negative and that 


can 
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In order to determine whether typical airfoils with neither sue- 




tion nor blowing have regions in which 


R 


cSRq 


> 0, that is, regions in 


which K >-0.925 (fig- 3).> the modified Schlichting method was used 
to calculate the distribution of K along “x for the MCA 64A010 and 
MCA 642AOI5 airfoil sections (reference 14) at zero angle of attack 

(figs. 8(a) and 9(a)). The distribution of K along x for v w =0 
at -all values of x Is independent of R c . For the MCA 64A010 airfoil 

section only the boundary layer over the forward 2 ^ percent of the air- 

Ra 2 R« 

^ D c ^ y c 


foil has 


R 


e 


dR* 


> 0. For the MCA 64£A015 airfoil, however, 


He 


1 Rfl 

for the forward 22^ percent of the surface. _ The region of — > 0 


8 r £ 

He^ 

Re 


> 0 


is greater for the MCA 642AOI5 airfoil section than for the MCA 64A010 
airfoil section because the MCA 642A015 airfoil section has larger 


values of 


dU 

dx 


over its forward portion than the MCA 64A010 section. 


(Compare figs . 10 and 11. ) 


It is to be noted that figures 3 and 6 indicate that by a suffi- 
cient increase in boundary-layer thickness the boundary layer over the 
forward portions of both the MCA 64A010 and the MCA 64^015 airfoil 

sections can be made stable. This conclusion follows from the result" 

that > 0 for K > -1 (fig. 6) and the result that 

SRe 


dK R0 



> 0 


for K > -0.925 
therefore make 


(fig- 3). 

R e c 

- — equal 

R e 


A sufficient 
unity unless 


increase in Rg (in 


U% 

dx 


is very large. 


&) can 
For 
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example, if the value of 9 at x = 0.35 on the NACA 64-A010 airfoil 
at Rc = 10T is multiplied by 3.^-5^ the value of Rg is increased 

Rg 

from 1221 to ij208 and the value of — 2. i s increased from 0.239 to 1.0. 

R 0 


When ^ = 0, K = -1 
dx 

plate flow, an increase in 
unstable. 


and 

R 0 




00. In this case, the flat- 


can only make the boundary layer more 


Effect of blowing at the stagnation point .- The analysis has shown 
that in a region of falling pressure on an airfoil having the usual type 
of surface (v w = 0) the disturbing effect of roughness particles can be 
decreased and the stability of the boundary layer increased by an 
increase in boundary- layer thickness. This result suggests that an 
effort be made to find a method for increasing the boundary- layer thick- 
ness without causing transition. The second part of the analysis inves- 
tigates such a method, namely, blowing in the region of the stagnation 
point. 

The work of reference 15 (see also appendix B) has shown that even 
for large blowing velocities the velocity profile in a stagnation-point 

flow, a flow in which U = f— ] x, does not have an inflection point 

\dx/s 

for y > 0. In fact, as v w becomes very large and positive (blowing) 
the velocity profile approaches the sine profile: 


u 

U 



which has an inflection point at y = 0. The value of Rg for this 

limiting velocity profile is 228 (appendix B). The value of R$ c at 

the stagnation point is therefore always greater than or equal to 228. 
At the stagnation point, however, Rg = 0. Therefore, Rg c is there 

always greater than Rg with the result that blowing at and near the 

stagnation point can bfe used to thicken the boundary layer without 
making it immediately unstable. 
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It is shown in appendix C that the Von Karman momentum equation 
upon which the Schlichting method is based is valid for a stagnation- 
point flow even when there are large velocities through the surface. 
The Schlichting method can therefore he used to calculate the develop- 
ment of the boundary layer near the stagnation point even when there 
are large velocities through the surface. 

As a result, blowing in the region of the stagnation point was 
investigated for both the NACA 64A010 and the NACA 64^015 airfoil 

sections at a = 0° and R c = 10^. The independent parameter was the 
quantity k which was given the values 0.709 and 7.09; these values 
are, respectively, 10 and 100 times the value of k at the stagnation 
point - without blowing. Because 



the factors 10 and 100 mean that the value of 8 at - the stagnation 

point is, respectively, \Jl0 and 10 times the value of 8 for zero 
blowing. The values of k-j_ and K are obtained from figure 12 and 

the chosen value of k. The magnitude of (v w ) g is obtained from 
figure 12 and the definition 



The behavior of the boundary layer is then computed by the modified 
Schlichting method (appendix A). 

At R c = 10^, a = 0°, and v w = 0 everywhere, the Reynolds num- 
ber Rq becomes larger than the critical Reynolds number R@ at -- 

0 

x = 0.038 for the NACA 64A010 airfoil section (fig. 8(b)) and at 
x = 0.060 for the NACA 64^015 airfoil section (fig. 9(h)). The 
boundary layer on both airfoils is thus unstable for the greater part 
of the region of falling pressure. When v w is made positive (blowing) 
near the stagnation point, Rg is increased and Rq c is decreased so 

that R 0 becomes equal to Rg c nearer the stagnation point than when 
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v w = 0. The values of Rg and Rg c , however, can then he made equal 
to one another for some distance by decreasing v w along the surface 
in the proper manner. When v w reaches zero, control of the boun dar y 
layer in the present examples ends. Of course the boundary layer can 
be made neutrally stable to the trailing edge by allowing v w to be 
negative (suction) from the point where v w = 0 to the trailing edge 
(see reference 6). In this case the boundary-layer thickness is first 
increased by the blowing and then decreased by the suction. The net 
effect depends on the value of x, the Reynolds number, the airfoil 
velocity distribution, and the magnitude of the blowing velocity at 
the stagnation point. 

The distribution of v w necessary to make Rg = Rg c ' is calculated 

by using the relation between Rg and Z together with the require- 
ment that Rg = Rg c . The relation between Rg and Z is 

Rg = U|/Zi/R^ 

but 

R e = R e c 

therefore 

Rg c = Ul/Z^; (18) 

In order to obtain the distribution of v w along x, the integration 
process described in appendix A is used. This procedure is briefly as 
follows: At x + Ax 


Z 


( 0 ) 

x+Ax 


= Z x + 



Ax 


The value of Rg^^ at x + Ax is then found from equation (18). The 

value of k(°2x is fouM by use of figure 1. The value of k^+Ax is 

found from Z^®) and hy use of the definition 

^ V ax /x+Ax 
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k = 



when and - 

k ( °) 

are known, 

x+Ax 

x+Ax 


and the functions 

S(K) 

and f(K). 


ki(0) is found from equation (4) 
x+Ax 

( 0 ) 


The value of 


I'dzy 
^dx/x+Ax 


is now 


found from equations (A2) and (A3). The iteration process described in 
appendix A is then used to find the final value of all the quantities 
at x + Ax. The value of L ~ v w at x + Ax is found from the equation 


f 


1 


^1 

/Z 


where 


f l = 

The procedure is repeated to find v w at x + 2&x and so forth. 

When v w becomes zero, the computation in the present examples is 
carried on with kp = 0 . 

The results of the computations for the NACA 61A010 airfoil are 
shown in figure 8 . The behavior of the curves in figure 8 (a) confirms 
the prediction of figure 6 that for v w = 0 an increase in "S 
increases K if K > -1. The distributions of v w along x for 
(v w ) 8 = 0.0151 and for (v w ) = 0.0515 are shown in figure 8 (c) for 

the region where _ v w > 0. The variation of_ Rg and Rq c for 
(v w ) s = 0 , (v w ) s = °* 01 51 ; and (v w ) s = 0.0515 is shown in fig- 
ures 8 (b), 8 (d), and 8 (e), respectively. In figure 8 (f) is shown the 
distribution of Re c /Rg along x for the three distributions of v w 

with x. The ratio of the boundary-layer thicknesses with blowing to 
those without blowing is shown in figure 8 (g). 

There are two main results of the computations for the NACA 64A010 
airfoil section. The first is that except for the region between 
x = 0.038 and x = 0.06 (fig. 8 (f)) the stability is decreased by 
the blowing; for x > 0. 06 the decrease is small. The second result; - 
is that, although there is a large increase in boundary- layer thickness 
near the stagnation point (fig. 8 (g)), the amount of the increase in 
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thickness decreases rapidly until at the pressure minimum the boundary- 
layer thickness for (v w ) g = 0.0151 is only about 2 percent and the 

boundary- layer thickness for (v w ) g = 0.0515 only about 5 percent 
greater than for (v w ) = 0. Thus, the only effect making transition 

less likely is the increased boundary-layer thickness and this increased 
thickness is present over only the foremost portion of the airfoil. 

The results of the computations for the MCA 64^015 airfoil are 

shown in figure 9- The behavior of the curves in figure 9(a) confirms 
the prediction of figure 6. The variation of v w along x for 
(v v ) 8 = 0.0111 and for (v w ) g = 0.0380 are shown in figure 9(c) for 

the region where vw >0. Figures 9(t>)> 9(<i)> and 9(e) contain the 
variation of Rg and of Rq c with x for (v w ) g = 0, (v w ) s = 0.0111, 

and (v w ) g = 0.0380, respectively. In figure 9(f) is Bhown the distri- 
bution of Rg c /Rg along x for the three distributions of (v w ) s 

along x. Figure 9(g) contains the ratios of boundary- layer thicknesses 
with blowing to those without blowing. There are again two main results. 
For the MCA 64-2A015 airfoil the ratio Rq c /Rq is greater for (v w ) g > 0 

than for (v w ) s = 0 for values of x between about 0.065 and 0.30 
(fig. 9 (f)). The ratio of boundary- layer thickness with blowing to 
that without blowing (fig. 9(g)) shows the same behavior as for the 
MCA 64-A010 airfoil. 

The chosen distribution of blowing has thus shown a greater effect 
in reducing the likelihood of transition on the MCA 6U2AOI5 airfoil 
than on the MCA 64A010 airfoil. On the MCA 6igA015 airfoil the chosen 
distribution of blowing has, in addition to the increase in boundary- 
layer thickness, also produced an increase in stability that is notice- 
able to about x = 0.30. 

Power required for blowing .- The power required for blowing is 


r QAH 

Eb 

where Q is the volume of fluid blown out per second, AH is the total 
pressure added to the fluid by a system within the aircraft, and Eg is 
the efficiency of the entire boundary-layer control system. 

The quantity Q AH is obtained by integrating v w AH around the 
entire airfoil surface] thus 
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Expressed as a drag coefficient, the power required for boundary-layer 
control is : 


-Dp = 


PEp 


— U q 3 c 
2 0 


d I 

E B 1 U 0 q Q c 


where Ep is the efficiency of the aircraft propulsive, system. Then, 


c % = r/ v »^ ax <19) 


If no fluid originates or is retained in the aircraft, the total drag 
coefficient is 



where 2 / 

\Cy 


is the wake drag coefficient' of the airfoil (reference lo) 
The expression for Cp^ (equation ( 19 )) can also be written as 




( 20 ) 



where 


m 


is an average pressure loss defined by 
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Therefore 



% /£h\ c 

E B\qo/ m 


Q 


The values of Cq for the NACA 6^-AOlO and NACA 64^015 airfoils are, 
for a = 0° and R c = 10^, 


Airfoil 

( y w) s 

C Q 

NACA 61*A010 

0.0151 

.0515 

0.00026 
. OOO58 

NACA 642A015 

.0111 

.0380 

. 00032 
. 00070 


Because of the small values of C 
appreciable only when the factor 


q, the drag coefficient Cp p will be 
Ep/£h\ 

— is a multiple of unity. 

E-D \ r» / 
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DISCUSSION 


The present analysis uses two approximate methods, namely, 
Schlichting’s method (reference 11) and Lin's approximate formula 
(reference 12). The tests of the Schlichting method in reference 11 
and the test' in appendix B for blowing, near the stagnation point give 
no reason to doubt that the Schlichting method correctly predicts the 
response of the boundary layer to changes in the variables that affect 
its behavior, namely, pressure gradient, Reynolds number, blowing, and 
so forth. The accuracy of the prediction, however, is known to depend 
on the particular case (see reference 11). 

Reference 17 investigated the approximate Lin formula and concluded 
that the predicted critical Reynolds numbers" agree well with those 
predicted by more elaborate calculations. Although the Lin formula 
has good accuracy, the exact prediction of the— critical Reynolds number 
of a velocity profile requires a precise knowledge of the first and sec- 
ond derivatives of the velocity profile u =~u(y). It therefore seems 
likely that the main uncertainty in the predictions of Rq^Rq by the 

combined Schlichting and Lin methods lies in the Schlichting method which 
assumes that all velocity profiles form the particular single -parameter 
family given by equation (Al), appendix A. 

In spite of the inexactness of the methods of analysis the result 
that an increase in boundary-layer thickness can increase the stability 
of-”the laminar boundary layer under the proper conditions is believed 
valid. It therefore follows that thinning the boundary layer will not~ 
always increase the boundary-layer stability." The present analysis, 
however, does not predict the effects obtained by placing slots on a 
surface (references 7, 8, and 9) because this investigation treats the 
effect of a change in boundary -layer thickness alone and slots change 
not only the boundary -layer thickness but also, by their sink effect 
(reference 9)> change the pressure distribution in their. vicinity and, 
by Introducing stagnation points, can produce convex velocity profiles 
on their downstream sides. 

It is noted that~the pressure drop along the surface required to 
make any thickening of the boundary layer result in an increase in sta- 
bility is not very large. For example, the NACA 61^015 airfoil at 

a = 0° has a large enough pressure drop over the forward 22^- percent 

of the surface to make any increase in boundary -layer thickness increase 
the stability. 

The result that the effective height of surface irregularities can 
be decreased by an increase in boundary -layer thickness without a decrease 
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in stability and the fact that the consistent attainment of laminar flow 
in flight is prevented by surface irregularities makes it important to 
find a method that can increase the boundary-layer thickness without 
itself causing transition. 


In the present work an investigation was made of the feasibility 
of increasing boundary- layer thicknesses by blowing over the foremost 
portion of an airfoil. In order to make the effects clear, large 
blowing velocities at the stagnation point were chosen. In order for 
the Von Karman momentum equation to be valid, however, it is necessary 


that v w be of the order of 


1 


or that 


• - (a- 


That is. 


either the boundary- layer-theory assumptions must be satisfied or (see 
appendix C) the flow must be a stagnation-point flow. For definiteness, 

-7==; at R c = 10 T , 

/rT 

- (S) .«• 


it is assumed that v w must be less than 


v w 


must 


thus be less than 0. 003 if U 


For the NACA 64A010 airfoil U = 17^-. 7x to .about x = 0. 003 with 
less than 10-percent error; thus the stagnation-point flow extends to 
about x = 0.003. In this region large values of v w do not invalidate 

the Von Karman momentum equation; therefore the Schlichtlng method 
remains valid for large v w for x less than about 0.003. For 

(v w ) b = 0.0151 however, v w > 0.003 for values of x between 0.003 
and 0.011 (fig. 8 (c)) and for (v w ) e = 0.0515, v w > 0.003 .for values 

of x between 0.003 and 0.012 (fig. 8 (c)). Therefore in the region 
lying roughly between x = 0.003 and x = 0.012 and about 0.01 in 
length, the values of v w are large enough to introduce an additional 

uncertainty into the predictions of the calculations for the NACA 
6UA010 airfoil. A similar region exists between x = 0.007 and 
x = 0.015 for (v w ) s = 0.0111 (fig. 9 (c)) and between x = 0.007 and 

x = 0.020 for (v w ) s = 0.0380 (fig. 9(c)) on the NACA 69^015 airfoil. 

The region in which the Schlichting method is probably invalid grows 
with an increase in (v w ) g so that it is inadvisable to make calcula- 
tions for larger values of (v w ) s 011 the NACA 64A010 and NACA 6kpA0 15 
airfoils at a = 0° and R c = 10?. 

Although the blowing velocities were so large that the boundary- 
layer assumptions were probably violated in a small area at the end of 
the stagnation-point-flow region, the boundary- layer thickness was 
increased significantly over only the foremost portion of the airfoils. 



32 


NACA TN 2752 


Whether larger blowing velocities can appreciably increase the boundary- 
layer thickness over larger portions of the "surface than found in the 
present work must therefore be determined experimentally. 


CONCLUDING REMARKS 


In order to investigate - the possibility that an increase in 
boundary-layer thickness could increase the local critical Reynolds 
number more than the local boundary -layer Reynolds number, computations 
have been made by combining the Schlichting method with the Lin method 
for the calculation of the critical Reynolds number of a velocity pro- 
file. The computations indicate that in a region of falling pressure 
on an impervious surface an increase in boundary-layer thickness can 
cause the velocity profile shape to be changed enough by the increase 
in effective pressure gradient so that the ratio of the local critical 
Reynolds number to the local boundary -layer Reynolds number is increased. 
It thus appears that~the local stability of the boundary layer can be 
increased by a local increase in boundary-layer thickness. The compu- 
tations also indicate that similar effects occur when there is flow 
through the surface; in this case the results depend on the effective 
flow through the surface as well as on the effective pressure gradient. 

These calculations suggest that an increase in boundary-layer thick- 
ness can decrease the disturbing effect of roughness particles without a 
decrease in stability. This conclusion is implied by the result that an 
increase in boundary-layer thickness reduces the velocity at a fixed dis- 
tance from the surface more than the change in velocity profile increases 
the velocity. One method of increasing the boundary-layer thickness, 
namely, blowing near the stagnation point, has been investigated theo- 
retically and seems to have limited potentialities because the blowing 
produces a significant increase in boundary-layer thickness only over 
the foremost portion of the airfoil. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., March 19, 1952" 



5A 


NACA TN 2752 


33 


APPENDIX A 

MODIFICATION OF THE SCHLICHTING METHOD 


The Schlichting method assumes that all velocity profiles of the 
laminar boundary layer form the single-parameter family 


§ = FjU) + KF 2 (tj) 


7 = 


5-, (x) 


(Al) 


where K is the velocity-profile shape parameter, 
and. F 2 are (see reference 11 ) : 


The functions 


F-l = 1 - e~^ 


and 


F 2 = Fi - sin ^ (0 < T) ^ 3 ) 


F 2 = -e 


-7 


(7 >3) 


The f un ction F^_ represents the asymptotic suction profile (reference 11 ) 
and the term sin ^ T|J represents an approximation to the Blasius pro- 
file for the flat plate (K = -l). With this relation for the velocity 
profile, Schlichting uses the momentum equation in the form 


Uf + f 2 + 5* \ dD _ T w 6 

V dx \ 6 ) v as V T7 tT 
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By introducing the substitutions 


k = £_ dU = z dIJ 
V dx bx 

kq = - ^ = f q \Jl 
V 


and 


T 9 
1 V 


pU 


where 


Z = 0 2 R 


cj 


the momentum equation can be written as 



The terms f and d*/e are functions of the parameter K which is a 
function of k and kq. “ _ _ 

In order to -compute the properties of the boundary layer, equa- 
tion (A2) is integrated step by step to find Z(x); therefore, G(k,kq) 
must be found at each step. Schlichtlng gives a plot of G(k,kq)j the 
computation of Z(x) therefore is made by using equation (A2), the 
definitions of k and kq, and the plot of" G(k,kq) (fig. 6 of 
reference ll). When any property of the shape of— the velocity profiles 
is needed, the distribution of the shape parameter K(x) must also be 
found. The parameter K, however, is given by Schlichting not as a 
function of k and kq but as a function of X and Xq; thus 
(equation (lb) of .reference II ) 
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\ + - 1 
1 - m( i - 1) 


(A^) 


where 


_ 2 

S 1 dU 

v dx 


and. 


, _ v wSl 

X 1 “ 

v 


In order to find K(x), it is therefore first necessary to find 
k(x) and kq(x) from k(x) and kq(x), the quantities available from 

the computation of Z. The computed distributions k(x) and kq(x) 
are converted to k(x) and kq(x) by two sets of curves, one, values 
of 1 against k for constant values of kq, the other, values of kq 
against k for constant values of kq. The distribution of the shape 

parameter K(x) is then found from a plot of equation (A^) (fig. 3 of 
reference 11 ). 


In the present investigation, properties of the boundary layer on 
an airfoil with blowing ^v w > 0^ had to be calculated. Because the 

curves given in reference 11 for use in calculations are confined to 
Vy % 0, it was necessary to extend them to v w >0. In doing this it 
was found that the curves of k(k,kq) and kq(k,kq) could be elimi- 
nated and the extension of the Schlichting method to v w > 0, simplified. 
The quantities k and kq then do not appear in a computation; K is 

found directly from k and kq. 


The quantities k 
of reference 11 for k 
obtain 


and kq are eliminated by solving equation (22) 
and equation (23) of reference 11 for kq to 
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and 



(AS) 


where g = g(K) =-»£-. When equations (A5) and (A6) are used with equa- 

5 1 

tion (a 4) the result is 


g 2 (K + 1) - g 


1 + 


l 1 - i> 


- k = 0 


(AT) 


Equation (20) of reference 11, however, is 


g 


1 + 


i 1 - §) K 


= f(K) 


Therefore, equation (A7) hecomes: 


g2(K +1) - fk 2 - k = 0 


(A8) 


a relation "between K, k, and k]4 The function K(k,kq), shown in 
figures 4(a) and 4(h), was calculated by fixing kq and calculating k 

for -2.099 = 0 from equation (A8). The functions g(K) and f(K) 

are given in reference 11. 


The curves of G(k,kq) given by Schlichting (fig. 6 of reference 11 ) 
can thus be extended to negative values of kq (that is v w > 0) by the 

use of equation (A3), figure 4, and the values of f(K), g(K), and 

b*/d(K) given in table 2 of reference 11. In the present work, how- 
ever, because of the large range of k and kq, it was more convenient 
to compute G(k,k x ) at every step in x by using figures 4(a) and 4(b) 

c* 

and curves of f(K), g(K), and — (K). 

8 


In order to Integrate equation (A2), the value of Z at x + Ax 
was found by the formula: 


(0) 

Z x+Ax 


= 



AX 
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_ ( 0 ) 

The value of Zx+Ax was used to find 
called Zx+Ax was "then found from: 


/dz\(°) 

\dx/ x+ Ax 


A new value of Z 


x+Ax 



This process was repeated either until there was no change between 

(if 


successive values of 


az\ 

^vx+Ax 


or to 


dZ \ 

^ /x+Ax 


When 


was 


not equal to 
was equal to 
was equal to 


'dZ 


A (l) 
V^/x+Ax 
'dZ\(l) 
y,dx/ x+ Ax 

’dZ\(°) 

/x+Ax 


x+Ax 

dZ\(2) 


the step length was reduced. When . . 

V^Vx+Ax 

the step length was not changed. When 

\ dx /x+ Ax 


the step length was increased. In no case was 


a value of Ax larger than 0.05 used. 


In order to begin a computation at the stagnation point, it is 
necessary to find the values of k and kq there. The requirement 

that ^ remain finite at the stagnation point (U = 0) means that 
G(k,kq) = 0 (equation A2). The values of k and kq at the stagna- 
tion point therefore satisfy the equation 


f 



(A9) 


The variables f and S*/S are functions of K only. Equation (A9) 
contains three variables, k, kq, and K. In order to find another 

relation between k, kq, and K, equation (A8) is used. There are then 
two equations In the two unknowns k and kqi 



kq = 0 
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and 


g 2 (K +1) - k - fkq = 0 


The solution is 


and 


- f 2 - g 2 ( K + 1) 

( 2 + fy - ^ 


(AlO) 



(All) 


Equations (AlO) and (All) give the values of k and kq at the stagna- 
tion point as functions of K. Because K at the stagnation point is 
usually unknown, it is better to have k and kq as functions of 
known quantities. Schlichting (reference 11) introduces a quantity Co 
defined as: 


C Q = 



By writing kq as 


*1 - - r % 

Rc \]j2 dx) 


and k as 
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is obtained. The quantities (v w ) s , Rqj and are usually known. 

The connection between k s , ki s , K s , and C Q is shown in figure 12. 

Although equation (A2) leads to 

dZ _ 0 ^ 

dx 0 


at the stagnation point, an application of L' Hospital's Rule results in 


a definite expression for 


(§Z\ . 

^dX/g’ 


it is 



z A ao + dfi ^ ao_ 

dx^ Sk dx ^1 



f 1 dG 

2 \fz Ski 


s 


For all computations in the present work ( — ] = 0 because for a 

Vdx J s 


symmetrical airfoil at zero angle of attack 

(**l\ 

the present work — — = 0. 

Vote L 


( d 2 U \ 
Vdxws 


= 0 and because in 
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APPENDIX B 

FLOW AT THE STAGNATION POINT FOR LARGE BLOWING QUANTITIES 

At the stagnation point the requirement that — remains finite 

dx 

although U = 0 means that G(k,k]_) =0 at the stagnation point. 

The expression for G(k,k^_) is 

G(k,kl) = 2jf - k(p + - k^j 

Therefore, at the stagnation point, 

f - *( 2 + t ) - k i - 0 

Now divide by k^; the result is 



For the stagnation-point flow, k is positive. The maximum 
value that f can ever have is 0.5; the maximum value of — is k.6b. 

e 

Therefore, when k]_ becomes very large and negative, equation (B3) 
becomes 


(Bl) 


(B2) 



(Bl*) 


The equation relating K, k, and k]_ also must apply; it is 


g 2 (K +1) - k - fk x = 0 


(B5) 
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Now divide "by kq; the result is 


g^K + 1) 

*1 



f = 0 


the maximum value of g^(K + l) is 0.25; therefore, when 
very large and negative, equation (B6) becomes 


-fc- = -f 

kl 


Therefore, for G = 0, k becomes large and positive as 
large and negative. If equations (b 4) and (B7) are combined, 
is 



Now (reference ll) 


f = jjL + ^1 - |-)k ( c q + cqK + C2K 2 ) 


and (reference ll) 


a* 


1 - ( g - -y 

C 0 + C]K + CgK^ 


Equation (B8) is then 


jl + ^1 - gK 2 (c Q + C]_K + c^g + 1 
By substituting the values (reference ll) 


-(a- 


K 


C1 - c 2 = 2 - - 


(B 6) 


becomes 


(B7) 


kq becomes 
the result 


(B8) 


= 0 (B9) 
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it“can "be shown hy trial that K = -1 is a root“of equation (B9) • 
The other two roots are found to he 

K = -4.95 " 

and 

K = 8.96 


Because both are outside the permissible range of K (-2.099 ^ K <£ 0) , 

K = -1 is the only valid root of equation (b 8) . Therefore, for a 
stagnation-point flow with very large positive v w the velocity profile 

becomes the Schlichting flat-plate profile 

S' sin (H <B10 > 

The relation between K, k, and k]_ at the stagnation point 
for all values of k^ is given by equations (AlO) and (All) and is 

shown in figures 4 and 12. By use of equations (AlO) and (All), it 
can be shown that k— ) 0, k^_— > 0.5, and K-> 0 as C Q 00 and that 

k — > 00 , kj — > - 00 , and K — ) -1 as C 0 — ) -«. From these results and 

figures 4 and 12 it is clear that for a stagnation-point flow -1 < K ^ 0. 

Therefore, the value of Rq c for a stagnation-point flow is always 

greater than 228 (see fig. l) and the velocity profile never has an 

inflection point for y > Oj velocity profiles for K > -1 have no 

inflection point for y > 0. 

The velocity profile predicted by the Schlichting method for a 
stagnation-point flow with large positive v w , namely, equation (BIO) 
is now shown to be an exact solution of the Navier-Stokes equations. 

The work of reference 15 has shown that the exact solution of the 
Navier-Stokes equations for a stagnation-point flow with large positive 
v w is 


u 

U 


sin 


\v w dxy 


The variable T) 


is 
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From the definitions of k and it follows that 

= _ v w 
k o dU 

® to 

In the limit for large v w (large k^_) , however, (see equation (b4) ) 



Therefore 


n = 



Now, from equation (17) of reference 11, 

-2- = c 0 + CjK + C2 K 2 
°1 

and from equation (19) of reference 11 

1 - (2 - |)k 

6 Cq + C]K + C2K 2 



Therefore 


a |cr\ 
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In the limit as v w gets large and posi-tive, K = -1; therefore, 



but, from reference 11, 


c o 



and 


c l - e 2 = 2 


6 

it 


Therefore, 


6 

5 1 



1 - 4 + 2 - 

it 


tl + 2 


6 = 6 

it « 


and in the limit 



Thus equation (BIO) becomes 


u 

U 



The Schlichting method- therefore predicts the correct velocity 

f n r u^ e sta S nation -P° int for large blowing. Note, however, 

hat the Schlichting method is incorrect in predicting that the velocity 
profile is the same as the flat-plate profile without blowing. 

V? t 5 e veloclt y Profiles predicted by the approximate 
Schlichting method (reference 11) with some exact-solutions of the 
boundary- layer equations for large v w in a stagnation-point flow 
(reference 18) is shown in figure 13 . The exact-solutions of the 

b ^? a T layer , e<3Uati ° nS f0r a stagnation-point flow are also solutions 
of the Navier-Stokes equations (reference 19 , p. 82 ). 
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APPENDIX C 


' / 

VON KAEMAN MOMENTUM EQUATION FOR STAGNATION-POINT FLOW 


The stagnation-point flow, the flow in the region in which 
U = ( — ) x, can he described by the Navier-Stokes equations in 

v<w s 

nondimens ional form; the equations are 


u Su + v^ii = -^ + — ^ 

Sy Sx R c \^ x 2 ^ y 2 j 


(Cl) 


u 


Sv 

Sx 


+ 



+ i/afv + s 2 v\ 

^ R c\Sx 2 Sy / 


(C2) 


The stream function for an inviscid stagnation-point flow with 
v ^ 0 at y = 0 is 

t = axy - bx 

( see reference 19 for the case b = 0) where 



tJ 0 c 


a 



and 


b 


_b_ 

Uo 



k6 
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The velocity components are U = ax and V" = -ay + b. By Bernoulli's 
theorem, the static pressure is 


P = P s 




For a viscous fluid the expression for the stream function is generalized 
to 


t = x0(y) - bx 

and the expression for the static pressure, to 

,2 


P = Ps - tt[x 2 + F(y,b)] 


The velocity components are 


Sty B$ 


v =--^=-0 + b 

Bx 


(C3) 


(Ch) 


The velocity components at y = 0 are u = 0 and v = b; at y = <» 
the x velocity component is the same as that for the inviscid fluid; 
that is 


u = U = ax 


(C5) 


It will now be shown that the Von Karman momentum equation results if 
equation (Cl) is integrated, with respect to y. From equations (c4) 
and the fact that the use of equations (C4) leads to a solution of-the 
stagnation-point flow with blowing (references 15 and 18), it follows 
that 


Bfu 

Bx 2 


= o 
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Equation (Cl) thus becomes 


Sx By Bx R c ^y2 


(C 6) 


Row integrate equation (C 6) with respect to y; thus, 


u dy + f v^dy = 


Bx 


dy 


0 


^ dy + i- F dy (C7) 
ox R C J 0 By 2 


where Z is a constant such that, for y > Z, u = U 
By integration by parts it follows that 

, ■I'- 

-l 0 

But 


= ax. 


r 1 du . 1 2 

I V — dy = vu I 
0 ay 


/ Bv j 

I u — dy 
0 ay 


^y 


Bv 


T= Jo ¥ dy + T " 


and, by continuity, 


therefore, 


Bv _ Bu 
By Bx 


£ v | ay = - u jf | dy + Dv » + ^ “i 4 * 


Equation (C7) becomes 

r z au 2 


Bu 


^ dy - U / g dy + Uv w = - 


Bx 


'0 




Jo 


(c8) 
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From equation (C3), 

dp 


dp q dp 

__ = -a^x; — is thus independent of y. By 

dx dx 


use of U. = ax, — can also he written as 
dx 


dp: = .y® 
dx dx 


(C9) 


Equation (C8) , after equation (C9) has been substituted and 
added and subtracted, is 






U 





dy + 



where 



for y > Z. After collecting terms and using 


(CIO) 


1 /du\ _ T w 

R c W/w pu 2 
K o o 

equation (CIO) becomes 



But 



= 9 
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k-9 


and 



(l - fry - S* 


therefore, the result is the Von Karman momentum equation 


dS 

dx 



v w T w 



Near the stagnation point (where 



✓ X 

the Von Karman 


momentum equation is therefore obtainable directly from the Navier-Stokes 
equations of motion and is thus valid there even for large v w . The 


Schlichting method is therefore not invalidated by the presence of 
large v w in a stagnation-point flow. 
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TABIE I.- SOEFLBMEKEAEJ IRFarattllOH FOR CURVES OF FIGURES 2 AUD 5 


(a] Values of J and related Information 


_ Sign of root In 

Value of J Range of kj. Range of K equation (6) 


Figures 2(a) end 5 (a), v w > 0 


10 * 

- 10 * < kx - -38 

-2.068 g K< 0 

10 

- 10.48 < ii < -1.071 

- 1.461 < K< 0 

1 

-L 366 g g -0. 1869 

- 1.111 < 1^0 

.1 

-0.2791 < g - 0 . 02113 

- 1.013 g Kg 0 

.001 

-0. 02287 g g -0 . 0002145 

- 1. 00013 g K g 0 

-.001 

-0.0122 g Iq. < 0 

-i. 924 gKg-i 

-.1 

-0. 1207 gljgO 

- 2 . 0 lBg K < -1 

-1 

-0. 38 < k x g b 

-2.099 g K g -1 

-10 

-1.20 g kjL < 0 

-2. 099 g K g -1 

-10* 

-38 < ki < 0 

- 2.099 g Kg -1 

10* 

-38 g ^ < 0 

- 2.068 g K g -1 

10 

-I.071 < k x g 0 

-1. 46 i g K g -1 

1 

- 0.1869 g Is^ g 0 

-Lin g Kg -1 


Figures 2(1>)' and 5(t>)j < 0 


1 

0.001 

0 S k i 5 0.02187 

-1 g K g 0 

2 

.1 

OJtii 0.1791 

-1 g K g 0 

3 

1 

0 < k-L < 0.3661 

-1 < K < 0 

4 

10 

0 <k x < 0.4770 

-1 g K g 0 

5 

10 * 

0 < i L < 0. 5000 

-1 g K g 0 

6 

- 10 * 

0 g k x < 0. 5000 

-1 g K < 0 

7 

-10 

0 g k x g 0. 3280 

-1< K< 0 

8 

-1 

0 < k x < 0.2639 

-1 < K < -0 9120 

9 

-1 

0.26395 k x g 0.4632 

-1.346 g K< -0.0120 

10 

-.1 

0.0218 < k x g 0.1277 

-L 785 5 Kg - 0 . 983 9 _ 

11 

-.001 

0.0002147 g k-L g 0.0123 

-L926 < Kg -0.9999 
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TABLE I. - SUPPLEHEJfCARY IflFORHATIGH FOR CURVES CF FIGURES 2 ASD ? - Concluded 
(b) Information for extreme values of K 



^Condition A: 1 + ^ 

f2j 


il/SK 


0; condition B: kj 
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Figure 1.- Variation of critical "boundary- layer Reynolds n umb er R with 

w c 

velocity-profile shape parameter K. 
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(a) T w > 0. 


Figure 2.- Variation of boundary- layer-stability parameter — w 

v w p 

u R<9 

velocity-profile shape parameter K for fixed values of boundary 
layer- control parameter J (see table I for values of j). The 
dashed line is the limiting curve (see equation (9)). 










(to) < 0. 

Figure 2.- Concluded 
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(b) 0.1 = k = 10. 

Figure 4-.- Concluded. 













ki e 0.111 


j | ■ | j Limiting value of K for k-]_ 
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Limiting value of K for k^ - 0 . 998 U 
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Figure 7.- Variation with velocity-profile shape parameter K of — _£. ; parameter of rate of change 

% dH e 

of roughness Reynolds number with boundary- layer thickness^ for fixed values of parameter of 

flow through surface k^. 




ON 
-F " 



(a) Variation of velocity-profile shape parameter K with position x along surface for NACA 64A010 
airfoil section at zero angle of attack for three values of (v w ) a , ratio of velocity through sur- 
face at the stagnation point to free- stream velocity. 

Figure 8.- Calculations for NACA 64A010 airfoil. 
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(b) Variation of boundary- layer Reynolds number Rg and, critical bo u ndary- layer Reynolds number Rq 
vith position x along surface for the NACA 64A010 airoil section for no flow through surface 
(v w = 0), a chord Reynolds number of 10 7, and zero angle of attack. 


Figure 8.- Continued. 





(c) Variation of v v , ratio of velocity through surface to free-stream velocity, with position x 
along surface for the RACA 64A010 airfoil section at zero angle of attack, a chord Reynolds number 
of 10^, and two values of (vw) s , value of v w at stagnation point, For (v w ) = 0.0515, 

0.04 s 

v w dx = 0.000287 and v w = 0 at x = 0.04j for (v w ) B = 0 . 0151 , 

and v v = 0 at x =» 0.04. 


/ 
J n 


V w dx = 0.000129 



Figure 8 .- Continued. 
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(a) Variation of boundary- layer Reynolds number Eg and critical boundary- layer Reynolds number 

°c 

with position x along surface for the NACA 64A010 airfoil section for (v w ) s , ratio of velocity 

through surface at stagnation point to free-stream velocity, equal to O.Olpi, a chord Reynolds 
number of 10^, and zero angle of attack. 

Figure 8.- Continued. 
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= 0.0515 | 

— |— v w variable 


(e) Variation of boundary- layer Reynolds number Rq and critical boundary- layer Reynolds number Rg^_ 
■with position x along surface for NACA 64A010 airfoil section for (v w ) s , ratio of velocity 
through, surface at stagnation point to free-stream velocity, equal to u.Opipj a chord Reynolds 
number of 10^, and zero angle of attack. 


Figure 8.- Continued. 
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(f) Variation of the ratio of critical bou n dary-layer Reynolds number Rq c to boundary-layer Reynolds 
number Eg with position x along the surface for the NACA 61t-A010 airfoil section at a chord 

Reynolds number of 10?, zero angle of attack, and fixed values of (v w ) g , the ratio of velocity 
through the surface at the stagnation point to the free -stream velocity. 


Figure 8.- Continued 
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(g) Variation with position x along surface of ratio of boundary -layer thic kne ss with blowing to 

O R 0)v v X) 

thickness without blowing — ~ for NACA 64A010 airfoil section at a chord Reynolds number 

Y 

of 10', zero angle of attack, and fixed values of (v v ) g , ratio of flow velocity through surface 
at stagnation point to free-stream velocity. 


Figure 8.- Concluded. 
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(a) Variation of velocity-profile shape parameter K with position x along surface for NACA 64-2A015 
airfoil section at zero angle of attack for three values of (v w ) s , ratio of velocity through sur- 
face at stagnation point to free-stream velocity. 


Figure 9.- Calculations for NACA 61^015 airfoil. 




_ayer Reynolds number % and critical boundary -layer Reynolds number R q q 
surface for the NACA 64-2A015 airfoil section for no flov through surface 

.ds number of 10^, and zero angle of attack. 

Figure 9.- Continued. 
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(c) Variation of v w , ratio of velocity through surface to free-stream velocity, with position x 
along surface for NACA 64gA015 airfoil section at zero angle of attack, a chord Reynolds number 

of 10 T, and two values of (v w ) B , value of v v at stagnation point. For (v w ) g = 0.0380, 


n 0.095 

nO.075 

/ v w dx = 0.000354 and v v = 0 at x = 0 . 095 ; for (v w ) = 0.0111, 

J 0 > 

/ v w dx = 0.000157 

J 0 

t r _ r\ v _r\ A7R 

CUU.12. v-y — \J a Li .A. — U,U| J m 


Figure 9.- Continued. 

; i 

1 






sag 


0.0111 


▼« variable 


(d) Variation of boundary-layer Reynolds number R Q and critical boundary-layer Reynolds number R fl 

“c 

with position x along surface for NACA 6^A015 airfoil section for (v v ) , the ratio of velocity 
through the surface , at the stagnation point to the free -stream Velocity, equal to 0.0111, a chord 
Reynolds number of 10?, and zero angle of attack. 


Figure 9 -- Continued. 
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(e) Vaxiutiuii uf boimdary-layer Reynolds number and critical boundary-layer Reynolds number 

” c 

with position x along surface for NACA 642A015 airfoil section for (v w ) s , ratio of velocity 
through surface at stagnation point to free-stream velocity, equal to O.0380, a chord Reynolds 
number of 10?, and zero angle of attack. ' 


Figure 9*- Continued. 
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(f) Variation of ratio of critical boundary-layer Reynolds number Rg^ to boundary-layer Reynolds 

number R g with position x along surface for NACA 6 ^A 015 airfoil section at a choird Reynolds 

number of 10^, zero angle of attack, and fixed values of (v w ) , ratio of velocity through the 

s 

surface at stagnation point to free-stream velocity. 


Figure Q-- Continued. 
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(g) Variation with position x along surface of ratio of boundary-layer thickness with blowing to 
thickness without blowing — — for NACA 6%3A015 airfoil section at a chord Reynolds number 

_ Wv-o 

of 10 f zero angle of attack, and fixed values of (v w ) s , ratio of velocity through surface at 
stagnation point to free-stream velocity. 


Figure 9 .- Concluded. 





Figure 12.- Variation at stagnation point— of velocity-profile shape 
parameter K ; pressure-gradient- parameter k, and parameter of flow 
through surface kp, with stagnation-point parameter C 0 . 










cJC 


8 : 



(a) No flow through the surface. C Q = 0. 


Figure 13.- Compar.' 
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Figure 13.- Continued. 
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